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Age and growth analysis is essential to fisheries management. Indirect methods to calculate growth are 
widely used; however, length frequency data analysis in sea cucumbers is complicated by high data 
variability caused by body wall elasticity. Here we calculated Isostichopus badionotus parameters of the von 
Bertalanffy growth function. In order to address bias produced by body wall elasticity, we compared the 
performance of four measurements and one compound index that combines different biometric 
parameters: the square root of the length-width product (SLW). Results showed that variability in length 
data due to body wall elasticity was controlled by using body length (Le) from the SLW compound index. 
Growth in I. badionotus follows a negative allometric tendency. Slow or zero growth periods were observed 
during October and November, when weather conditions were adverse. 

Most of the targeted holothurian species by fishers through the world are reported as over-exploited 1,2 . A 
recent worldwide analysis reveals that 20% of fisheries were found to be depleted, 38% over-exploited, 
14% fully exploited (i.e. no potential margin for expansion) and only 27% moderately or under- 
exploited 3 . In order to maintain stocks at levels capable of producing maximum sustainable yield, countries 
should adopt appropriate measures based on available scientific data 4 . Holothurians need particular attention due 
to their vulnerability to overfishing, biological characteristics and commercial value 1 . 

In recent times some species are being targeted in Latin America, where scientific knowledge is limited. 
Specifically, ten years ago, a large Isostichopus badionotus 5 stock was identified on the northern Yucatan 
Peninsula continental shelf 6 . This stock has attracted the interest of sea cucumber traders and fishers in devel- 
oping this potential new fishery. However, lack of regulations and scientific knowledge on the species in the region 
has hindered development of any fishery sustainability guidelines, threatening the stock with overexploitation. In 
the western Atlantic, I. badionotus has been captured in Cuba 7 , Venezuela and Panama 8 . To date, studies on this 
species have focused on distribution, abundance and biomass assessment 6,9 , however, other biological character- 
istics, such as age and growth have not yet been studied. 

Isostichopus badionotus is widely distributed in the Atlantic, from Central America to West Africa, and from 
Bermuda to Ascension Island 10 . This species is found at depths of 0 to 65 m on shallow muddy-sand and sea grass 
beds 11 , bare hard surfaces, coral rubble, dead coral, pavements, erect macroalgae and turf algae 12 . Few data on 
biometric relationships are available for I. badionotus 6,s,UM . 

Most methods address growth through length frequency data analysis, which depends on body measurements 
and biometric relationships. Precision measurements of weight and length in sea cucumbers are challenging. 
Body wall elasticity, digestive system content and coelomic fluid produce bias in size frequency data analysis 15 . 
This challenge can be addressed by using three alternative set of variables: 1 ) weight, immersed weight 16,17 , drained 
weight 15,18,19 and/or gutted weight 13,15,20 ; 2) length, either contracted 21,22 or relaxed 13,23 ; and 3) compound indices 
that combine different biometric parameters 17,24 30 ; the latter one allows calculation of basal area and volume that 
can then be applied to generate more precise biometric relationships. 

The objective of the present study was to document biometry and growth for a I. badionotus populations of the 
northern Yucatan Peninsula continental shelf (Mexico). Dorsal length, total weight and gutted weight were 
measured. A compound index and its transformation were analyzed and compared, as were the resulting 
biometric relationships. Isostichopus badionotus growth parameters were calculated by indirect methods. 
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Results 

Morphological measurements and population structure. A total of 
2347 sea cucumbers were measured during the study period, and five 
frequency distributions analyzed from the body length (L), width 
(W), total weight (TW), muscle wet weight or gutted weight (GW) 
measurements, and the compound index SL W (the square root of the 
length-width product) (Fig. 1). Mean body length ± SD (25.3 ± 
5.2 cm; Fig. la) was almost twice mean body width (14.2 ± 
3.1 cm; Fig. lb), with sizes ranging from 3 to 45 cm and a 
population mode of 23 cm. The SLW frequency distribution 
showed the minimum coefficient of variance (18.6%; Fig. lc), 
whereas the maximum coefficient of variance (56.9%) was 
observed in the GW frequency distribution (Fig. le). The SLW 
distribution was multimodal and ranged from 2.4 to 28.3 cm. Both 
the TW and GW frequency distributions were unimodal (mode = 
600 g and 400 g, respectively). The TW frequency distribution had a 
mean of 552 ± 263 g with weights ranging from 1 to 1800 g(Fig. Id). 
Weight loss from TW to GW was 37.5% (Fig. If) and, contrary to 
expectations, coefficient of variation (CV) was higher for GW than 
for TW. 

Biometric relationships. The regression equation between L and 
SLW used to generate estimated length was: Le = 1.11 + 1.28 
SLW (r 2 = 0.75, p < 0.001). The overall Student t-test showed 
that I. badionotus grows allometrically (ft 3, p < 0.05), with 
negative allometric ft values ranging from 1.5 to 2.7 (Fig. 2). 
However, the length-weight relationships by location showed [> 
values near 3 (isometric value) at Progreso. In all cases, the 
relationship between body length and weight was significant (p < 
0.001). The Le measurement explained almost 80% of variation in 
body weight and was a more effective predictor of weight than L 
(Table 1). This is because Le was derived from two measurements 
(square root of length multiplied by width) rather than one, which 
improved the accuracy of body weight estimates (Figs 2c and 2d). 



The equation (1) produced poor body weight estimates at lengths 
< 18 cmand>30 cm when using length alone (Figs. 2a and 2b). The 
L-GW relationships showed a point cloud (points inside the dotted 
circle in the Figs. 2b and 2d) that does not fit the general tendency of 
most observations. This point cloud corresponded to measurements 
recorded in 2009, 52% of which were taken in September. 

Analysis of mean condition factor (Kn ± SD) showed that Kn 
values were lower in September 2009 than in other months, overall 
and by location. Mean Kn values generally ranged between 0.37 ± 
0.25 in September 2009 and 1.98 ± 0.25 in November 2010 (Table 2). 
However, when the data were broken down by location, mean Kn 
values were higher during September 2010. 

Growth analysis. The ELEFAN I method (Electronic Length Fre- 
quency Analysis) yielded the best estimated growth parameter results 
both in the pooled sample and for Celestun: Ie„ = 31.6 cm, k = 
0.6 yr -1 (overall); Le x = 31.9 cm, k = 0.5 yr' 1 (Celestun) (Table 3). 
For Sisal and Progreso, the best results were produced by the 
PROJMAT method (Projection Matrix): Le m = 38 cm, k = 
0.2 yr~' (Sisal); Le« = 23.5 cm, k = 0.7 yr" 1 (Progreso). 
Asymptotic gutted weight (GWJ by location was 599 g for the 
pooled sample, 638 g for Celestun, 929 g for Sisal and 192 g for 
Progreso. Parameter C was 1 for Celestun, 0.7 for Sisal, 0.6 for 
Progreso and 0.8 for the pooled sample. These results indicate 
intense intra- annual oscillations in growth. Winter point (WP) 
parameter values were 0.8 (overall and Celestun) and 0.9 (Sisal and 
Progreso), implying minimum growth was present in October 
(overall and Celestun) and November (Sisal and Progreso). 

The likelihood ratio tests indicated differences (~/ 2 test, df = 3, p < 
0.001) in growth between locations. All comparisons showed differ- 
ences, both in isolated VBGF parameters (y 2 test, df = 1, p < 0.015) 
and when combining Le x and k (x 2 test, df = 2, p < 0.015). The 
exception was the Sisal vs. Progreso comparison. In both the isolated 
Le„ and k parameters, as well as their combination, the values were 
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Figure 1 | Frequency distributions of measurements and a compound index for sea cucumber Isostichopus badionotus. A total of n = 2347 sea 
cucumbers I. badionotus were measured during the study period in 2009 and 2010 off the Yucatan Peninsula, Mexico, (a) L distribution; (b) W 
distribution; (c) SL W distribution; (d)TW distribution; (e) GW distribution; and (f) weight loss from TW to GW of n = 1165 sea cucumbers, loss is 
indicated as a percentage and shown in bold. L: dorsal body length; W: dorsal body width; SLW: square root of length multiplied by width; TW: total 
weight; GW: gutted weight; SD: standard deviation; and CV: coefficient of variation in percentage, lowest value is shown in bold. 
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Figure 2 | Relationships between estimated lengths and weights from sea cucumber Isostichopus badionotus. Points inside the dotted circle are data 
recorded during 2009. TW: total weight; GW: gutted weight; and Le: estimated length from the equation Le = 1.11 + 1.28 SL W, where SL Wis square root 
of length multiplied by width. 



not significant (p > 0.015). These two locations may have similar 
environmental features, while those at Celestun are different. 
However, the t a value showed a significant difference between Sisal 
and Progreso, whether in isolated Le^ and k or when combined. 

Graphic representations of growth (Figs. 3a and 3b) suggest that 
sea cucumbers younger than one year in age had Le < 18.9 cm (GW 
< 151 g) in Celestun, Le < 9.8 cm (GW < 30 g) in Sisal, and Le < 
15.3 cm (GW < 56 g) in Progreso. The largest recorded specimens 
were about five years of age in Celestun (Le > 30 cm and GW > 
400 g), ten years of age in Sisal (Le > 35 cm and GW > 700 g), and 
four years of age in Progreso (Le > 22 cm and GW > 140 g). These 
results indicate that individuals measured during the study period 
were older than two years at both Celestun (mean ± SD: 25 ± 4.1 cm 
and 304.2 ± 158 g) and Progreso (mean ± SD: 20 ± 2.9 cm and 
124.1 ± 81.1 g), and older than four years at Sisal (mean ± SD: 
26.3 ± 4.8 cm and 350.8 ± 10.2 g). 



Table 1 | Length-weight relationship parameters by location. Le: 
estimated length in cm; TW: total weight in grams; GW: gutted 
weight in grams; a: is a constant; ji: allometric coefficient; SE: 
standard error; r 2 : coefficient of determination; p: test signifi- 
cance level; and n: sample size 



Locality Relationship 



tSE) 



(i(±SE) 



Celestun Levs. TW 
ievs. GW 

Sisal ievs. TW 
ievs. GW 

Progreso ievs. TW 
Lews. GW 



0.18 ± 0.02 
0.04 ± 0.01 
0.1 1 ±0.02 
0.09 ± 0.02 
0.03 ± 0.01 
0.02 ± 0.03 



2.48 ± 0.04 
2.76 ± 0.08 
2.61 ±0.05 
2.54 ± 0.07 
2.94 ±0.14 
2.86 ±0.37 



0.77 «0.001 1350 
0.78 «0.001 591 
0.81 <0.001 882 
0.78 <0.001 460 
0.82 <0.001 115 
0.40«0.001 115 



Discussion 

Biometric relationships are useful to carry out transformations of 
length to weight or vice versa, and then calculate growth parameters. 
Length and weight relationships have not been established for all sea 
cucumber species, and this is the first study that applies a compound 
index on I. badionotus. Length, as well as width, is another important 
biometric characteristic. Body length is a poor estimator of body 
weight in sea cucumbers, and width provides little useful data per 
se. However, when length and width are combined to produce the 
SLW index, a modified version of the method proposed by Yamana 
and Hamano 25 , the analyses improved. 

Using estimated length (Le) from SLW enhanced weight estima- 
tion accuracy by 80% in L. badionotus. This improvement in body 
weight estimation might be due to a significant variance reduction in 
length and width during organism handling. When measurements 
are carried out on live sea cucumbers, the act of handling can cause 
them to change shape elastically. As a result, the organism can 
deform along different symmetry axes (e.g. from left to right, top 
to bottom and forward to backward), changing individual measure- 
ments. Other advantages of SLW are that it facilitates in situ sea 
cucumber length and width measurements, is less time consuming 
to take measures, and thus lowers organism stress level. The SLW 
compound index also showed multimodal frequency distribution, 
another advantage when identifying cohorts in growth studies 
through indirect methods. 

The growth parameters reported here for I badionotus are similar 
to those reported for other species in the family Stichopodidae 
(Table 4). A study done by Herrero-Perezrul et al. 19 in the Pacific 
Ocean (Mexico), reports a growth rate for Isostichopus fuscus minor 
to the one found for I. badionotus in our study. This difference could 
be due to various reasons: 1) the deep water temperature in the East 
Pacific is colder than in the Gulf of Mexico, 2) I. fuscus specimens 
were measured underwater and, according to our field experience, 
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the contraction of animals is inevitable, even underwater. 3) Herrero- 
Perezrul et al. 19 used a simple version of the VBGF, while we used a 
version that considers seasonality. Different authors have shown the 
existence of shrinking in sea cucumbers during certain periods as 
consequence of seasonal growth patterns due to changing envir- 
onmental conditions 17,31,32 . As above, it is important to consider the 
version of the von Bertalanffy growth function with seasonality for 
sea cucumbers. 

Isostichopus badionotus exhibits a negative allometric growth, and 
our results provide strong evidence that growth in this species vary 
seasonally. The Eq. (2) results showed that growth was slower during 
October and November (WP = 0.8 and 0.9, respectively). Seasonal 
growth variations in the pooled sample showed growth to increase 
about 80% (C = 0.8) in the peak "growth season" (June to 
September) with an immediate 20% decline in October. In conclu- 
sion, the I. badionotus population on the northern continental shelf 
of the Yucatan Peninsula displays a variable growth pattern, with a 



rapid growth rate prior to October and a slower growth rate during 
October and afterwards. 

Intra- annual oscillations in I. badionotus growth may be related to 
environment (e.g. water temperature, salinity, dissolved oxygen, and 
transparency), as well as to biological processes such as seasonal 
cycles in feeding behavior for energy storage 16 . Growth variations 
coincided with seasonal changes (e.g. rain and cold front seasons) 
on the Yucatan coast, so these fluctuations and food availability may 
in turn impact condition factor. Growth rate was enhanced when 
environmental conditions were favorable and the organisms accu- 
mulated energy for an eventual reproduction season. In the present 
study, relative condition factor (Kn) did not decrease during the slow 
or zero growth seasons; on the contrary, it improved, perhaps due to 
nutrient storage in the body wall. There are evidences that physico- 
chemical properties of the surface sediments linked to food availabil- 
ity for deposit-feeding detritivores determined the distribution 
pattern of Australostichopus mollis juveniles 33 . When organic matter 
availability was low, A. mollis juveniles ceased feeding and moved 
faster over the sea floor to new areas where the level of organic matter 
in the sediment was sufficiently high to stimulate feeding 34 . 

Another characteristic related to seasonal growth in I. badionotus 
are inter-annual changes in biological condition. Relative condition 
factor was higher in 2010 than in 2009, principally during September. 
Gutted weight (GW) gave important information about biological 
condition of organisms. GWhad a high coefficient of variation (CV 
= 56.9%), caused by individuals measured in 2009, as shown in the 
L-GW relationships (Fig. 2). The large variations in GW indicate that 
body walls were thinner in 2009 than in 2010; a trend confirmed by 
lower relative condition factor (Kn) values in 2009, particularly in 
July and September (Table 2). These variations occurred from June to 
November. Similar to condition factor, growth rate decreased to the 
point of body shrinkage during periods of environmental stress (e.g. 
algal blooms, hurricanes and winter storms), most likely in response 
to limited food availability. 

Hurricanes have been shown to have short-term impacts on mar- 
ine environments 35 . In September 2010, tropical storm Matthew 36 
impacted Belize and southern Mexico, producing torrential rainfall 
on the northern Yucatan coast in September and October. It was 
during September 2010 that the highest Kn values were recorded. 
Another phenomenon that may have affected sea cucumber stocks in 
this period was a large harmful algal bloom in summer 2008 off the 
Yucatan coast 37 . It negatively affected food availability in the regional 
marine ecosystem, and caused widespread mortality in different 
marine organisms groups. Indeed, sea cucumbers disappeared from 



Table 3 | Growth parameters of Isostichopus badionotus. St = best 
estimated obtained through either ELEFAN I or PROJMAT applied 
to all samples;^ = jackknifed estimate; A % = percent difference 
between Stand the average of all Jackknife pseudovalues; CVj = 
coefficient of variation of jackknifed estimate; Le^ = estimated 
asymptotic length (cm); k = growth rate (year - '), t a = time at 
which Le„ is zero (years); <p' = Comparative growth. Rn and 
SSQ = Goodness of fit indices for ELEFAN I and the sum of 
squared differences for PROJMAT, respectively 
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Figure 3 | Individual growth and relative ages by location of sea cucumber Isostichopus badionotus. Left: seasonalized VBGF (von Bertalanffy growth 
function) by location in cmyr~'; Right: seasonalized VBGF by location in gyr -1 . Relative ages were generated for each Leand GWdatum by replacing the 
growth parameters calculated by ELEFAN I and PROJMAT with their respective values in the seasonalized VBGF, using sub-routines programmed in 
Wolfram Mathematica 9®. Le. Estimated length; GW: gutted weight. 



a number of locations and abundance was decimated in the northern 
and western coastal areas 38 . Small I. badionotus individuals were 
collected in summer 2009, an indication of recruitment and stock 
recovery. The thinner body walls recorded in 2009 may therefore be 
the aftermath of the widespread mortality and food scarcity caused 
by the harmful algal bloom. 

The sea cucumber body wall is an important organ for nutrient 
storage during feeding periods 39 . A laboratory study showed that I. 
badionotus enters a 12-month dormant phase in response to adverse 
environmental conditions. Coelomic fluid protein concentration 
increased, suggesting optimization of the protein synthesis process 40 . 
A similar phenomenon has been reported for A. japonicus during 
dormancy periods, when body wall tissue crude protein, lipids and 
carbohydrates contents reach maximum values 41 . These authors 
found that body wall biochemical composition was significantly 
influenced by the environment, particularly by seasonal changes. 
This has been observed also in other species of sea cucumbers of 
the family Stichopodidae, like A. mollis, whose growth and survival 
rates fluctuate significantly and are suboptimal at poor-quality envir- 
onments or low food availability 42 . 

Our results confirm that I. badionotus has developed physiological 
mechanisms to face adverse conditions, for instance with a radically 
reduction of metabolism 40 . Nonetheless, our data also suggest that 
the species still remains sensitive to meteorological and oceano- 
graphic phenomena due probably to changes of water temperature 
and food availability. Zamora and Jeffs 43 have showed that sea 
cucumbers can respond negatively to an increase in seawater tem- 
perature in terms of feeding, metabolism and growth. 



Even though most sea cucumber fisheries are managed without 
age and growth studies, understanding these parameters can help in 
improving management measures. The SLW index reduces the vari- 
ability of body length, which facilitates the calculation of growth 
parameters by indirect methods. The seasonal growth patterns in 
sea cucumbers can be well approached by one seasonally oscillating 
version of the VBGF, which may have implications in fisheries man- 
agement. The data provided here will be vital to fishery managers, 
who need to consider environmental variability and biological stock 
condition when developing sustainable management strategies. 
Furthermore, these results are potentially relevant to develop aqua- 
culture systems for this and other similar sea cucumber species in 
tropical ecosystems. 

Methods 

Organism collection and measurement were done at three locations off the northwest 
coast of Yucatan state, Mexico: (i) 20° 51' 48"N, 90° 24' 00" W (off Celestun); (ii) 
21° 09' 15" N, 90° 02 09" W (off Sisal); and (iii) 21° 17 11 N, 89° 39' 07" W (off 
Progreso). Isostichopus badionotus individuals were collected during four months in 
2009 (July, September, November and December) and 1 1 months in 20 1 0 (all months 
but March). Unfavorable weather prevented collection in August and October 2009, 
March 2010, and May to December 2010 at Progreso. Individuals were collected at 10 
to 18 m depth by SCUBA diving. The smallest individuals (from 3 to 12 cm length) 
were taken between 15 and 18 m depth on rocky reef areas with encrusting algae. The 
largest individuals were found on rocks, pavements and sandy beds with turf algae 
and hard corals (mainly blushing star and rose corals). Once collected, individuals 
were relaxed by placement in shaded containers with 10 L seawater and 24 g mag- 
nesium chloride (MgCl 2 ) for 5 to 10 min. They were then measured. 

Length (L), width (W) and total weight (TW) were measured for each individual. 
Length was measured dorsally from the anus to the center of the tentacle crown 
(curved length). Width was also measured dorsally, at the widest point 17,29 , excluding 
papillae on the ventrolateral margin and the ventral surface. Both length and width 
were measured to the nearest 0.1 cm with a tailor's tape. Total weight was measured 
with a digital balance (Torrey LEQ-5/10 Scale, ±1 g). In some individuals, the gut, 
gonad (if present), respiratory trees and coelomic fluid were removed and body wall 
weight (muscle wet weight or gutted weight) (GW) measured. Body length and width 
were used to calculate the compound index SLW (square root of length multiplied by 
width) 25 . Coefficients of variation were calculated for this index and each body 
measurement parameter. Collections were granted by fishing license no. 
DGOPA.08263.290709.2467 from CONAPESCA (Comision Nacional de 
Acuacultura y Pesca). 

Biometric relationships. Lengths were recalculated (Le) using the regression 
between L vs. SLW to increase growth parameter accuracy 25 . Description of relative 
changes in morphology was done using the allometric equation: 

y = ax>' (1) 

Where a and /5 are parameters, x is length (L or Le), and y is weight (TW or GW). 
Furthermore, four simple regression analyses were applied: L vs. TW, L vs. GW, Le 
vs. XWand Le vs. GW. Parameters were calculated with non-linear regression 
analysis, an effective tool in allometric studies since it simultaneously considers 
variation in both x and y variables 44 . The loss functions were evaluated with 
Marquardt's compromise algorithm 45 , which minimizes the sum of squares (SSQ) of 
the difference between sample weights and pseudo-weights 46 . Convergence was 



Table 4 | Growth parameters of /. badionotus and others species 
in the family Stichopodidae. L, = asymptotic length (cm), k = 
growth rate (yr _1 ), q>' = Comparative growth. *Le x > = estimated 
asymptotic length (cm) 



Species 


L (cm) 




<P' 


Stichopus vastus 54 


31.58 


0.55 


2.74 


Stichopus chloronotus* 5 


34.20 


0.45 


2.72 


Isostichopus fuscus ]9 


36.12 


0.18 


2.37 


Astichopus japonicus 55 


36.70 


0.33 


2.65 


Thelenota ananas 56 


66.30 


0.20 


2.94 


Isostichopus badionotus (Pooled sample, 


31.60* 


0.50 


2.77 


present study) 








Isostichopus badionotus (Celestun, present 


31.90* 


0.60 


2.71 


study) 








Isostichopus badionotus (Sisal, present study) 


38.00* 


0.20 


2.46 


Isostichopus badionotus (Progreso, present 


23.50* 


0.70 


2.58 



study) 
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reached at 50 iterations in each regression. A Student t-test was conducted to test the 
hypothesis of isometry (/? — 3). Allometric growth occurred when either ft > 3 
(negative) or < 3 (positive). Additionally, organism condition was analyzed using 
the relative condition factor 47 Kn as follows: Kn = GW/GW, where GW is gutted 
weight of individual sea cucumbers and GW is predicted gutted weight. The Eq. (1) 
was applied to predict GW as follows: GW = a(Le) /i . 

Growth analysis. Different Le class intervals were evaluated and two indirect 
methods were applied: 1 ) Electronic Length Frequency Analysis (ELEFAN I) 48 ; and 2) 
Projection Matrix (PROJMAT) 49 . One seasonally oscillating version of the von 
Bertalanffy growth function ( VBGF) was applied as follows 50 : 

U =Loo 1 1 — e~ ~ ' 0> + © ski2n{l " ' l} ~ ® sin2Tt( '° ~ ' s) l } (2) 

Where L t is predicted size at age t; L rJD is asymptotic size; k is curvature parameter per 
year -1 , expressing the rate at which L rJO is approached; and t Q is theoretical "age" if 
organism were to have a size equal to zero. Considering seasonal growth oscillations, 
t s defines the start of the convex segment of a sinusoid oscillation with respect to t = 0, 
and C is the parameter reflecting the intensity of seasonal growth oscillation. To 
calculate Le^, Le was used rather than L rJ> . 

Growth parameters were generated with the Length Frequency Distribution 
Analysis (LFDA) Version 5 software package (http://www.fmsp.org.uk/Software. 
htm), and compared using the growth performance index 51 <p'. The method was 
chosen with lowest variation <p' between intervals, and the interval with the best 
goodness-of-fit. The influence of each monthly sample in estimating growth para- 
meters was determined with the estimator A% and a jackknifed coefficient of vari- 
ation CVj, as described by Levi et al. 52 . Based on the growth coefficient /?, estimated 
from the allometric equation GW — a (Lef, the Eq. (2) was also expressed in terms of 
GW 46 . 

Relative ages at certain lengths and weights were generated by replacing the growth 
parameters with their respective values in the Eq. (2), using sub-routines pro- 
grammed in Wolfram Mathematica 9® (license number 3422-6765). 

Growth for individuals from the three sampled locations was compared by like- 
lihood ratio tests 53 , and three simultaneous tests (Celestun vs. Sisal, Celestun vs. 
Progreso and Sisal vs. Progreso) were performed to evaluate pair-wise differences in 
growth. These tests were equivalent to testing h independent hypotheses at a sig- 
nificance level. By applying Bonferroni inequality to the three tests resulted in an 
experiment wise significance level < h m a. An a = 0.015 was chosen, resulting in a 
0.045 experiment wise significance level. First, a simultaneous comparison was done 
of all three parameters under the null hypothesis: Ho: L mJ = L^ 2 \ Ki — K 2 ; t ol = t o2 . 
The model selection process was then done by sequentially altering the number of 
parameters in the comparison. 
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